On 
ON 



IX, 



Interfacial adsorption in two— dimensional Potts models 

E. Car Ion 

Laboratoire de Physique des Materiaux, Universite Henri Poincare, Nancy I, 

B. P. 239, F-54506 Vandoeuvre-les- Nancy Cedex, France 

F. Igloi 

Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, Hungary 
Institute for Theoretical Physics, Szeged University H-6720 Szeged, Hungary 

W. Selke, F. Szalma 
Institut fur Theoretische Physik, Technische Hochschule 
j^jL " D-52056 Aachen, Germany 

The interfacial adsorption W at the first-order transition in two-dimensional q-state Potts 
models is studied. In particular, findings of Monte Carlo simulations and of density-matrix renor- 
malization group calculations, at q = 16, are consistent with the analytic result, obtained in the 

i-Q . Hamiltonian limit at large values of q, that W oc t~ ' on approach to the bulk critical temperature 

T c , t — \T C — T\/T c . In addition, the numerical findings allow to estimate corrections to scaling. 
Our study supports and specifies a previous conclusion by Bricmont and Lebowitz based on low- 

£Nj ■ temperature expansions. 
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Q\ ' 1. Introduction 

ON ■ The interface between two phases may become unstable against the appearance of a third phase. This wetting 

phenomenon has been studied in various circumstances, considering different materials and geometries, both experi- 
mentally and theoretically!!! 

Good candidates for modelling wetting are g-state Potts modelscrcl, where the two phases correspond to distinct 
boundary states, say, 1 and q, at opposite sides of the system. In that case one observes an excess adsorption of the 
q non-boundary states at the interface between T' rich and V rich domains. 

o ■ 

In particular, in two dimensions, the interfacial adsorption, W, is found to diverge on approach to the critical 
temperature T c like W oc t~ u , where t = \T C — T\/T c . For q— 3 aad 4, the bulk transition is continuous and u can be 
J^ ■ expressed in terms of canonical bulk critical exponents, u> = v— j3n For larger values of q, the transition is of first order. 
On general, phenomenological groundsa, one may expect then u> = 1/3, as observed in Monte Carlo (MC) simulations 
on the two-dimensional Blume-Capel model.El However, previous simulations on the two-dimensional 20-state Potts 
model suggested that uj may be significantly larger than 1/3, albeit a pronounced curvature in the corresponding 
log-log plots for W(t) was noticed indicating that the asymptotic behavior may not have been reached.cl A reason- 
able explanation of the numerical findings was then offered by Bricmont and Lebowitz.13 Based on low-temperature 
expansions, the critical region is argued to be extremely narrow in the Potts case, in contrast to the Blume-Capel 
case.B However, a numerical confirmation remained to be done. 

Motivated by a recent intrigiung analysis of the interfacial tension of the i«-state Potts model in two dimeHsions 
employing a field-theoretic representation for configurations of the interfaced (extending prior considerationsu) , we 
decided to reconsider the somewhat irritating problem on the value of u> for q > 4. Indeed, advances in methods and 
computer facilities allow now to explore numerically the critical region more deeply than some years ago. In addition, 
we dealt with the issue analytically by studying the Hamiltonian limit of the Potts model for large values of q exactly. 

The outline of the article is as follows. In the next Section, the model is defined, and the numerical methods of our 
choice, MC simulations and density-matrix renormalization group (DMRG) calculations, are introduced. Then, we 
discuss the results obtained from those methods. In Section 4, the analytical findings in the Hamiltonian limit are 
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presented. A summary concludes the paper. 

2. Model and numerical methods 

The two-dimensional g-state Potts model is described by the Hamiltonian 

summing over neighboring sites {ij),(i'j') on a L X M rectangular lattice, with the Potts variable ny = l,2,...q. 
In the thermodynamic limit, L, M — > oo, the model displays a bulk phase transition at the critical temperature 
ksTc/J = l/ln(y/g + 1)0, where ks is the Boltzmann constant. The transition is of continuous type at q < 4, while 
it is of first order at larger values of q. 

To introduce an interface, we add a column of M sites on the left side boundary, fixing the Potts variable to be in 
the state rii, = 1, and another column on the right side boundary with tir — q. The Potts variables on the top and 
at the bottom of the lattice may be connected by periodic boundary conditions ('periodic case'). Alternatively, the 
interface may be pinned at the lower and upper boundaries by inserting there additional boundary rows with fixed 
states being T' on the left half and V on the right half (pinned case). 

By examining typical Monte Carlo equilibrium configurations below T c , it is seen that an excess of non-boundary 
states is adsorbed at the interface, as illustrated in Fig. 1. This wetting phenomenon may be described quantitatively 
by the net interfacial adsorption W per unit length 

W = 1/M J2 (< <U,„ y >l: Q - < S nb , rHj >l:l) (2) 

(ij),nb 

summing over all L x M sites in the inner part of the system; the < > brackets refer to thermal averages; the index 
nb denotes the non-boundary states, nb = 2, ..., q — 1; the subscripts l:q and 1:1 refer to systems with corresponding 
fixed states at the boundaries, i.e. with and without interface. The net adsorption W is closely related to the density 
profiles 

n bi:b2 (s,i) = 1/M ^ < ^«,n« >6i:6a ( 3 ) 

i 

describing the variation of the density of state s by going from the left side, fixed in state b±, to the right side, fixed 
in state &2, of the lattice, summing over each column, j = 1, ...M, with i running from 1 to L. Obviously, 

W = ^2(n 1:q (nb,i)-n 1:1 {nb,i)) (4) 

i,nb 

We computed profiles apd net adsorption numerically, using Monte Carlo techniquesa and the density-matrix renor- 
malization group method.113 In the MC simulations, the standard single-variable flip algorithm was applied, for system 
sizes L < 256 and M < 2000, studying mostly the periodic case, augmented by a few runs for the pinned case. To 
obtain accurate equilibrium data, we typically averaged over several realizations (using different random numbers), 
each consisting of 10 6 Monte Carlo steps per site, including at least 10 5 steps for equilibration. Error bars resulted 
from averaging over the ensemble of realizations. Usually, we set q = 16_syhere the bulk phase transition is strongly 
first order with a bulk correlation length at T c of a few lattice spacings.Eil To search for possible g-dependences, we 
also simulated models with q= 20 and 40. Furthermore, because clusters, formed by neighboring sites in the same 
state, seem to play an interesting role in the wetting phenomenon, see Fig. 1, we also computed cluster properties in 
the MC study. 

In the DMRG method, one considers strip-like lattices, i.e. M — ► oo, while L is finite. The algorithm, introduced 
by White ia 1992113 for the study of the low-lying spectrum of quantum spin chains, has been extended in several 
directionsJij In the present study we follow Nishino'sEj formulation of the method adapted to treat classical two- 
dimensional systems, where the approach is used for constructing iteratively approximate transfer matrices, starting 
from strips of small width (say, L = 8) which can be also handled numerically exactly. At each DMRG step the strip 
width is increased, and the configurational space is truncated efficiently through the projection into smaller subspaces 
with the help of appropriate density matrices. 



We do not need to describe details of the widely used DMRG procedure here: a good introduction, together with 
recent developments, can be found in Ref-Jlfl|. The DMRG technique has already been applied to the two-dimensional 
Potts model, both at continuous, q < 4,0-ffij and first-order transitions,^ > 413. In the former case, bulk and surface 
critical exponents have been calculated with a high degree of accuracyuJ. In the present study we are interested in 
rather large values of q, where the standard DMRG method is less suitable. An appropriate, powerful variant, the 
"pseudo - spin" version, has been introduced in Rcf. [16|. It enables one to treat systems with special values of q, 
namely q = p 2 , p = 2, 3, . . .. We used that algorithm in the current work. Most calculations were performed at q = 16, 
augmented by some calculations at the q = 9, for strips of widths up to L = 80. We kept up to m = 60 states pet 
block with a typical truncation error of e~ 10" 6 . From the dominant eigenvector of the transfer matrix the profileslild 
and the interfacial adsorption W were calculated. 

3. Numerical results 

The crucial quantity, which we computed numerically, is the net adsorption W as a function of the system size, 
L, M, and temperature t, W(L, M, t). The main aim is to analyse its critical behavior as t — > or T — > T c , in the 
thermodynamic limit L, M — > oo. For that, one may extrapolate the numerical data for W to that limit, or one may 
study systems for which finite-size effects can be neglected. We also tried to perform a finite-size scaling analysis on 
W(L,M,t]B, see below. 

In any event, accurate data are needed. Their quality can be conveniently tested by comparing results obtained 
from the MC simulations and the DMRG calculations, as illustrated in Figs. 2 and 3. 

In Fig. 2, profiles, Eq. (3), are shown for the 16-state Potts model with and without interface, demonstrating again 
the excess adsorption of non-boundary states at the interface. In Fig. 3, the increase of the net adsorption W with 
increasing width L of the Potts model, at various temperatures, is displayed. In both figures, the finite-size effect 
arising from the length M of the system can be disregarded. In fact, M is infinite in the DMRG approach. In the simu- 
lations (periodic case), M was checked to be sufficiently large so that W approached closely W(L,M = oo, t), with the 
characteristic crossover value depending, of course, on the width L and the distance from criticality t = \T C — T\/T c . 

As exemplified in Figs. 2 and 3, data from both methods do, indeed, agree nicely, being obviously accurate and 
reliable. Slight systematic deviations seem to become significant only for quite wide systems, say, L > 60. A reason- 
able extrapolation of our data for W to the thermodynamic limit is feasible for temperatures T < 0.998T C , see Fig. 
3. Closer to T c , both numerical methods would require much larger system sizes demanding extremely large storage 
and/or computing time. 

In the thermodynamic limit, one expects W oc t~ w as t — ► 0. To monitor the approach to the asymptotic behavior, 
one may consider the effective exponent 

u e ff(t) = -dhiW/dlnt (5) 

with ijj = u> e ff(t = 0). In Fig. 4, numerical estimates of the effective exponent are shown. At t > 0.02, we determined 
U) e ff from MC data for systems being large enough to disregard finite-size effects. Actually, finite-size effects are much 
stronger in the pinned case than in the periodic case, and only MC data for the latter case are included, with L = 64 
and M > 400. At t < 0.01, estimates are based on extrapolating MC and DMRG results for W to the thermodynamic 
limit, see Fig. 3. Note that the net adsorption is calculated at discrete temperat ures, tj, and the effective exponent 
may be approximated by uj e ff(t) = — hi(W(ti)/W(ti+i))/ ln(ti/ti + i), with t= ^(tiU+i), and ii+i < U. Error bars 
in Fig. 4 stem from a proliferation of the error in the net adsorption W . Evidently, to e ff depends strongly on the 
distance from criticality t. For instance, in the range 0.1 > t > 0.002, it changes from about 1.2 to about 0.5 when 
moving towards T c . Accordingly, an average critical exponent in that interval would be supposedly significantly larger 
than the true asymptotic critical exponent u), as observed beforeE Presuming u) = 1/3, see Refs. 4 and 6 as well 
as Section 4, it follows from Fig. 4 that the asymptotic regime is very narrow, and corrections to scaling are quite 
pronounced. To quantify these corrections, one may postulate the standard ansatz 

W(t) = W t- 1 / 3 (l + at x + bt 2x + ...) (6) 

The coefficients can be calculated from chi-square fits to the numerical data near T c . Using, for instance, a fit to 
the points in the interval 0.02 < t < 0.1 (shown in Fig. 4), leads to a net adsorption W reproducing very well the 
numerical findings both closer to T c and further away from T c . Eventually, systematic deviations show up when 
further lowering the temperature, reflecting the need for additional correction terms in Eq. (6) in that region, see 
Fig. 4. The exponent x characterising the corrections to scaling is found to be rather small, x — 0.14 ±0.06. Because 



of its smallness, we included the leading, with the exponent x, and subleading, 2x, terms in the ansatz (6). The error 
bars arise from using a variety of plausible fitting intervals and points. 

From MC simulations of Potts models with larger number of states, q — 20 and 40, one may conclude that the 
corrections to scaling, at t > 0.01, are rather insensitive to the concrete value of q. 

From general considerationsQcl, one expects two diverging lengths at the interface, in the direction parallel to the 
interface, £||, and perpendicular to it, £±, with £j_ oc t~ u and £|| ex t~ 2u> , as t — ► 0, sea also Section 4. For strip-like 
systems, M — ► oo, the following finite-size scaling expressions can be then motivatedQ 

W{L,M = oc,t) = t- u> w 1 {W) (7) 

in the limit of Lt^ >> 1, and 

W(L,M = oc,t = 0)ocL (8) 

for L >> 1. Indeed, the DMRG results indicate that Eq. (8) seems to be satisfied rather well already for strips of 
moderate width, say, L > 16, see also Ref. 3. On the other hand, the numerical data do not suffice to establish the 
scaling form (7) with u> — 1/3. In fact, a 'reasonable' collapse of our data on an apparent scaling function w\ might 
be achieved with a somewhat larger value, uj k, 1/2. Careful inspection, however, reveals that systematic deviations 
from a unique scaling curve w\ set in for large arguments Lt^ . Thence, larger systems close to T c had to be studied 
(which are, at present, out of reach, because of limitations in storage and computing time). Indeed, u) w 1/2 is merely 
an upper bound of the true value of 10. This observation corroborates the above mentioned finding on Lo e ff'- One 
has to include corrections to scaling to demonstrate consistency of the numerical data with the theoretically expected 
asymptotic behavior, due to the narrowness of the critical region. 

In closing the Section on the MC and DMRG results, we remark that the thermally averaged largest cluster of 
non-boundary states, as computed in the simulations, seems to diverge by approaching T c from below. A detailed 
analysis would be desirable, but it is beyond the scope of the present study. 

4. The Hamiltonian limit for large values of q 

In the .following we consider the Hamiltonian limit of the Potts model with strong vertical and weak horizontal 
couplingsO. The transfer matrix in the vertical direction has the form T = exp(— H), with the one-dimensional 
Hamiltonian HtB 

L-l L <j-l 

i=l i=\ fc=l 

ni is the Potts variable on site i and M* denotes the flip operation M* \rii) = \i%i + fc,mod q). The strength of the 
transverse field, h, at the transition point is h c = l/q. Quantities of physical interest are derived from the ground 
state, l^o), and from the energies of the ground state and the first excited state, Eq and E\, of (9). 



The Hamiltonian limit of the Potts model has been treated recentlyt3 for free boundary conditions. In that case, 
the solution has a remarkably simple form in the vicinity of the transition point for large values of q. Repeating the 
same type of considerations for models with an interface, fixing the variables at the boundaries in the states T' and 
V (l:q) (as before), one finds that the ground state sector of the Hamiltonian (9) is spanned by the vectors 

\ipij) = 1 11 .. .ln b n b . ..n b qq...q) , (10) 

where a non-boundary state is given by \n b ) = \/ ^fq{\2) + |3) + . . . + \q — 1)), and the positions of the domain walls 
separating the boundary and non-boundary states are denoted by i(= 1, 2, . . . , L) and j = (i, i + 1, . . . , L). The 
diagonal matrix-elements of these states, (ipi,j\H\i/Jij) = —L — (j — i)t, with t(= hq — 1, \t\ <C 1) being the distance 
from the critical point, are smaller by an amount of O(l) compared to any other states, like those containing boundary 
states in the domain of non-boundary (nb) states. The Hamiltonian in the ground state sector, spanned by the vectors 
(|l0|), can be written as 

H g = -(j - i)t - h^q-[(a+ + a-) + (b+ + b~)] , (11) 

up to a constant; the operators, a ± and b , which move the positions of the domain walls in (J10|), are defined as 



a± \^i,j) = \ipi±i,j), l<i<j; a + \il) iti ) = a \ipij) 
b \4>i,j) = \i>i,j±i), i<j<L; b~\ipi,i) = b + \ip itL ) 








(12) 



In the continuum limit, when L> 1, i> 1 and j 3> 1, but x = i/L = 0(1) and y = j / L = 0(1), the Hamiltonian 
O) can be written in the form of a differential operator 



H g ip(x,y) 



hy/g f d 2 d 2 

L 2 \dx 2 dy 2 



+ tL(y - x) 



ip(x,y) = Eip(x,y) , 



with the boundary condition < x < y < 1. 

At the critical point, t = 0, the solution of the eigenvalue problem ( |13| ) reads 

tp{x,y) = 2[sin(7rfcia;) sin(7rfc2?/) — sin(7rfc2x) sm.(wkiy)] , 



(13) 



(14) 



with fci = 1,2,... and fc2 = fci + 1, fci + 2, . . .. For the ground state, ipo(x, y), one has k\ = 1 and &2 = 2, whereas for 
the first excited state fci = 1 and fc2 = 3. Thus the energy gap is AE = Ei — Eq = 5ir 2 h^/qL~ 2 , and the correlation 
length parallel to the interface behaves as £|| ~ (AE)^ 1 ~ i 2 . Since the correlation length perpendicular to the 
interface is, at the critical point, limited by the width of the system, £x ~ L, one arrives at ^|| ~ Q_, in agreement 
with the form mentioned above. 

The density profiles satisfy the relations ni :q (l,x) — n\-. q (q, 1 — x) and n\ :q {x) — ^Jii: q {nb,x) = 1 — ni-. q (l,x) — 

nb 

n\ :q (q,x). From the ground state, one obtains 



ni-. q (l,x)= / dx' / dy[V'o(x,y)] 



■K 1 . 

— ( 1 — x) H — sin 2nx 

2 4 



— ( 1 — x ) H — sin 4irx 
2 y ' 8 



sin itx 

9 



(15) 



For small x, one finds Tii :g (l,a;) = 1 — (107r 2 /3)x 3 + 0(x 6 ), whereas for x close to one, the profile behaves like 
ni :q {l, x) ~ (1 — x) 10 . The profile of non-boundary states is symmetric and its maximal value is given by ni :g (l/2) = 
1/2 + 32/(9tt 2 ) = 0.86025. 

For the (1 : 1) boundary condition, the profile of non-boundary states tends to zero for large values of q. Therefore 
the interfacial adsorption W , Eq. (2), at the transition point may be approximated by 



W f MA l 

— = J ni: q {x)dx = - 



35 

72tt 2 



= 0.3826 . 



(16) 



Thus, at the transition point, W is, indeed, proportional to L, see Eq. (8). Note that the prefactor seems to depend 
on q, being, at q = 16, about 0.3, according to the DMRG calculations. 

Below the critical point, t ^ 0, we consider the eigenvalue equation (O) in terms of the new variables x+ = (y+x)/v2 
and X- — (y — x)j\[2. Then 



hy/q2 ( d 2 d 2 



L 2 



\dx\ 



dx 2 



tLV2x- 



ijj{x + , X-) — Exp(x + , X-) 



(17) 



with the boundary condition < x_ < x + < \j\[2 and < x_ < \/2 — x+ < l/\/2. Now the eigenfunction ip can be 



written as ij){x-\-,xJ) — (f>-\-{x+)4>-{x-). 4>+(x + ) satisfies the free-particle equation —d 2 4>+/dx 2 i _ 
(f)^(xJ) is the solution of the Schrodinger equation of a particle in a linear potential 



L 2 E+/(2h^q)^ 



h^/q2 d 2 
L 2 dxl 



-tLV2x- 



f>_(x_) = E_4>{x-) 



(18) 



Equation (Hg) leads to bound states, and the energy scale (both for the ground state and the excited states) is set 
by e ~ t 2 l. Hence the temperature dependence of the .parallel correlation length is given by £|| ~ (AE)~ 1 ~ t~ 2 / 3 , 
in accordance with the phenomenological considerations!! On the other hand, the extent of the bound states sets the 
length scale £j_ ~ t^ 1 ^ 3 , which is then proportional to the interfacial adsorption, W. Thence, in the Hamiltonian limit 
for large values of q, one has w = 1/3. 



5. Summary 

In this article, critical interfacial properties of two-dimensional g-state Potts models at the bulk first-order phase 
transition have been studied. We applied two numerical methods, Monte Carlo simulations and the density-matrix 
renormalization group approach, mainly at q — 16. Furthermore, we considered analytically the model in its Hamil- 
tonian limit at large values of q. 

The different methods lead to a consistent description of the critical behavior. The interfacial adsorption W diverges 
on approach to the phase transition temperature as W ex t~ w , uj = 1/3, with a very narrow asymptotic region. The 
strong corrections to scaling are characterised by a small exponent, x = 0.14 ± 0.06, which seems to depend (if at all) 
only weakly on the number of Potts states, q. At the critical point, W diverges linearly with the width of the system 
(being indefinitely long in the direction parallel to the interface). The proportionality factor has been calculated in 
the Hamiltonian limit. 

The value of the critical exponent u, u) = 1/3, is typical for wetting phenomena at bulk transitions of first or- 
der in two dimensions. It follows from general, phcnomcnological considerations, based on an effective interface 
HamiltonianQ, as well as from calculations on various microscopic multi-state models, such as Potts and the Blume- 
Capel models. The strong corrections to scaling and the narrowness of the,-critical region are, on the other hand, 
features which are specific for Potts models. They have been predicted beforeQ by using low-temperature arguments, 
and they have been quantified in this study. 
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FIG. 1. Typical Monte Carlo equilibrium configuration of the two-dimensional 16-state Potts model with an interface, at 
T — 0.99T C . The '1' domain is on the left hand side, and the '16' domain on the right hand side. Shown are clusters of distinct 
states. A MC system of size L — 60 and M = 60, periodic case, was simulated, but only a part is depicted. 



FIG. 2. Profiles, summed over the non-boundary (nb) states, ni-.b(i) = y~]ni:b (nb,i) for systems with interface, b = 16, and 

nb 

without interface, b — 1, at T = 0.99T C and L = 24, as obtained from the DMRG method (open circles) and MC simulations 
(periodic case with M — 400; full diamonds). 

FIG. 3. Net adsorption, W, vs. inverse width, 1/L, of the lattice at (from bottom to top) t = 0.01, 0.005, 0.002, 0.001, and 
0.0005, depicting DMRG (open circles) and MC (periodic case with M < 1000; full diamonds) data. 

FIG. 4. Temperature, t, dependence of the effective exponent of the net adsorption W , u> e ff(t), for numerical data 'free of 
finite-size effects', see text. The dashed curve corresponds to the fit to Eq. (6), with Wo = 3.280, a = —1.977, b = 0.939, and 
x = 0.165, quantifying the corrections to scaling. 
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